Massless Dirac fermions, gauge fields, and underdoped cuprates 
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We study 2+1 dimensional massless Dirac fermions and bosons coupled to a U(l) gauge field 
as a model for underdoped cuprates. We find that the uniform susceptibility and the specific heat 
coefficient are logarithmically enhanced (compared to linear-in-T behavior) due to the fluctuation 
of transverse gauge field which is the only massless mode at finite boson density. We analyze 
existing data, and find good agreement in the spin gap phase. Within our picture, the drop of the 
susceptibility below the superconducting Tc arises from the suppression of gauge fiuctuations. 
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Recent experiments have indicated the existence in the 
normal state underdoped cuprate superconductor of a 
gap with the same anisotropy as the d-wave supercon- 
ducting gap. One proposed explanation involves spin- 
charge separation: an electron in these highly correlated 
materials is a composite object made of a spin 2 neutral 
fermion (spinon) and a spinless charged boson (holon). 
The suppression of normal state magnetic excitation seen 
in NMR and neutron scattering is thus viewed as a sin- 
glet pairing of neutral fermions in the absence of coher- 
ence among holons. As a possible realization of this idea, 
two of us have taken the t-J model (which is believed 
to capture essential physics of Cu02 planes) and devel- 
oped a slave boson mean field theory |l| that extends the 
local SU(2) symmetry at half-filling to the finite concen- 
tration of holes by introducing a SU(2) doublet of slave 
boson field. Among the mean field phases reported in 
Ref. ^ the so-called staggered flux (sF) phase (which 
is connected to d-wave pairing phase by a local SU(2) 
transformation) was argued to describe the pseudogap 
in underdoped cuprates. The low energy physics of this 
phase can be described by massless Dirac fermions, non- 
relativistic bosons, and a massless U(l) gauge field which 
together with two massive gauge fields forms SU(2) gauge 
fields that represent the fluctuations around the mean 
field. 

The purpose of this paper is to address the low en- 
ergy effective theory of the sF phase as a U(l) gauge 
theory problem. Although Dirac fermions coupled to a 
gauge field had been considered in several contexts in 
the past y], we shall see that interesting new physics 
emerges when massless Dirac fermions are coupled to a 
gauge field that is also coupled to a compressible boson 
current. More specifically, the Lorentz symmetry break- 
ing due to coupling to the bosons results in the renor- 
malization of fermion velocity which have consequences 
on physical properties such as uniform susceptibility Xu 
and electronic specific heat c^ . Experimentally, Xu of 
underdoped cuprates begins to decrease with lowering 
of temperature far above the superconducting Tc P-pl . 
Electronic specific heat experiments [p|J7| show that 7(r) 



(= cf{T)/T) behaves quite similar to Xu- Although con- 
stant Wilson ratio (7/x«) is a hallmark of Fermi liquid 
theory, the anomalous temperature dependence calls for 
a departure from the time-honored theory of most met- 
als. We make a case that the puzzling normal state be- 
havior of Xu and 7 may be viewed as enhancement over 
linear-in-T Xu and 7 of Dirac fermions due to logarithmic 
decrease of Dirac velocity caused by fermion-gauge field 
interaction. 

We begin with the following continuum effective La- 
grangian for our problem 

C = *as(9^7'' + iaM7^)*a. + 

b* {do -tiB+iao)b- -^b* (V + imfb. (1) 

2mB 

The Fermi field *„« is a 2 x 1 spinor: *|^ = (/j^^^, f^^^), 
v|/t^ = (/|^o,/2se)i where a = 1,2 labels the two Fermi 
points, s — l,..,N labels fermion species {N — 2 for 
physical case s =t,i), and e,o stands for even and odd 
sites, respectively. The 7'^ matrices are Pauli matrices 
(7°,7\7^) = {<^^,cr\a^) and satisfy {7^,7''} = 2S'"' 
{fi,iy = 0,1,2). *c<s = ^L7°- In the sF phase of 
Ref. |l|], the fermion dispersion near the fermi points is 
anisotropic, but we rescale it to an isotropic spectrum 
_E(k) = f£)|k| where vd = \JvfVi, the geometric mean of 
the two velocities (v2 is proportional to the energy gap) . 
We set ud = 1, unless otherwise specified. The gauge 
field a^ = (flQ, a) corresponds to the a?' part of the SU(2) 
gauge fields of Ref. |]1|. The terms in Eq. (|l|) involving 
the Bose field b (representing charge degree of freedom) 
are believed to play several important roles, including 
the suppression of chiral symmetry breaking (Neel order- 
ing PI) and instanton effects |q|. Most importantly, the 
compressible boson current screens the oq field, making 
it massive. Unfortunately we do not have a detailed un- 
derstanding of our boson subsystem. Therefore we shall 
draw upon only a few of qualitative features of the Bose 
sector while focusing mainly on the Fermi sector of the 
theory. 

Eq. (n^) carries certain similarity to the uniform 
resonating valence bond (uRVB) gauge theory fPIJTC 



proposed to describe optimally and slighly overdoped 
cuprates, and some of the theoretical framework can 
be carried over to our problem. As in the uRVB case, 
the internal gauge field a^ does not have dynamics of 
its own, but it acquires dynamics from the polariza- 
tion of fermions and bosons. Integrating out the matter 
fields generates the self energy term for the gauge field 
^a = 5Q'A'(n^'' + I\!^)av, up to quadratic order. The 

given by 



fermion polarization W^ from the two Dirac points is 



nri?) 



kn 



(2^ 



ir[GF{k)rGF{k + qh''], (2) 



where Gpik) = -~{ik^'y^)~^ is the fermion Green's func- 
tion and fc, q denote 3-momentum; for example, k — 
(fco = (2n -I- l)7rr, k). In the Coulomb gauge, the spatial 
part and the time part of the gauge field are decoupled, 
the propagators being D°^{q) = (n°9(<7) -|-n°;"(<7))-i and 
D''{<1) = i5ij-qrq^/q^)D^{q) {i,j = 1, 2), with i?^(q) - 
{Ilp{q) + W-j^{q))~^. As mentioned earher, the bosons 
should have a finite compressibility (11^(9 ^- 0) 7^ 0) 
so the time component of the gauge field becomes mas- 
sive (at finite temperature I]Pp{q — > 0) is also nonzero 
and contributes to the screening of ao field), but the spa- 
tial part of the gauge field, which is purely transverse, 
remains massless even at finite boson density and tem- 
perature, as long as the bosons are uncondensed (as in 
the spin gap phase). In the remainder of this paper, we 
will focus on the effect of this massless mode, ignoring 
the ao field. 

In the absence of detailed understanding of the Bose 
sector, we assume that the transverse gauge propaga- 
tor is dominated by the fermion part. In other words, 
D^{q) « Dp{q) = 1/Iip{q). This approximation, which 
is often used in the uRVB gauge theory, may not be fully 
justified in our case, but it allows us to organize the 
infrared behavior of our theory within 1/A^ expansion. 
The full expression for analytically continued transverse 
polarization function I{p{Lo,c\) at finite temperature is 
rather complicated and therefore we shall not write it 
here, although it is used later in the evaluation of gauge 
fluctuation contribution to Xu and d^ . In the limiting 
case of T > |q| > \lo\, we have 

ni;(^,q)«-zCi^ + C2^, (3) 

while in the zero temperature limit, 

Imn^(u;,q) = ~Ns\gn{Lo)e[\u\-H)^uj^ - c^ /^ 
ReIlj^{u,q)^Ne{\ci\~\u:\)^^^Z^/8. (4) 

To the leading order in 1/N , fermion self energy due 
to transverse gauge fluctuations is 

^(^) ^ ^ E / (§^^^^^(^ + ^h'^F (9), (5) 



wheie D^p{q) = [Sjj — qiq^ / c^) /H-^ . At zero temperature, 
the self energy is ||ll|] 

/■ ifq ka + go 



8 ^ ' W (2^ 



)Hk + q)^V^' 



. ^ [ (fiq {k^ + qx){ql-ql)-'2qxqy{ky + qy) 
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(27r)3 



q2(fc + qY\/q 



yl dPq {ky + qy){ql-ql) -2q^qy{k^+q^) 



(27r)3 q2(^ + ^)2^ 

We find, for |k| > |fco|, 

S(fc) = -cikofAf){k) + 2cik- ^Ai{k) (7) 

with c = 4/(37V7r2) and Ao{k) w ^i(fc) w ln(A/|k|), 
where A is a UV cutoff. Now the pole of the renormalized 
Green's function Gp{k) = (Gpik)^^ — S(fc))^^ occurs at 



ii;(k) = |k|(l-4/(A^7r2)ln(A/|k|)). 



(8) 



Note that the presence of compressible bosons and the re- 
sulting breaking of Lorentz symmetry is crucial to have 
logarithmic velocity renormalization. Indeed, in the ab- 
sense of bosons, the gauge propagator (gauge indepen- 
dent part) is given by D'^'^iq) = 8/N{S^^~q^q^/q^)/ ^f 
and the zero temperature fermion self energy takes the 
form S — ik^j'^f{k'^); therefore the velocity is not renor- 
malized. 

Treating the quasiparticles described by Eq. (||) as 
"free", we calculate Cy and Xu up to 0{1/N°): 

c^' = (9/7r)C(3)A^T'(l + (S/A^tt^) ln(A/T) + ..) 

Xu = (2/7r) ln(2)A^ T{1 + {8/Nt:^) ln(A/r) + ..), (9) 

(C(3) = 1.202). These results are believed to be valid for 
two reasons: 1) Im^o,i('^ + iO+,k) = for |i^| < |k|, so 
the quasiparticles are well-defined. 2) Unlike the usual 
Fermi liquid theory, the free particle response function 
vanishes as T — > 0. To the extent the Landau parame- 
ters in Fermi liquid theory enter as in mean field theory, 
this means that Landau parameter correction vanish in 
T — > |I3]. Indeed, it will be shown shortly that the 
calculation of Xu and c^' from the free energy shift due 
to gauge fluctuation yields the same results. 

The enhancement of c^' seen here finds its counterpart 
in the more familiar problems such as electron-phonon 
interaction in metals |14|, uRVB gauge theory pffl, and 
half- filled Landau level |13], where interactions induce 
mass enhancement which manifests iteself in the specific 
heat. In the nonrelativistic analogues, however, mass 
renormalization does not necessarily result in the en- 
hancement of compressibility and uniform susceptibility 



|4,n6 17 1, because the corrections are tied to the Fermi 



surface [Q. The crucial difference in our case is that 
there are only Fermi "points" instead of Fermi "surface" . 
Thus in contrast to the nonrelativistic case, we find that 
the susceptibility is also renormalized such that the Wil- 
son ratio "f{T)/xu{T) is constant. In fact, the Wilson 
ratio is the same as that of free Dirac fermions because 
quasiparticles are well-defined and Fermi-liquid type cor- 
rections are absent, as discussed earlier. 

To check this conclusion, we calculate Xu and c^' in 
a gauge invariant way, using the correction to the free 
energy due to gauge fluctuations. We consider only the 
leading correction in \/N , which is 0(1/7V°): 
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AF~- — — d^q dioniuj) tan 
{2ttYJ J_^ 



_i/ Imn^(co,q) 

Ueni.(t^,q) 



(10) 



The entropy shift AS (= ~dAF/dT) due to gauge fluctu- 
ation has two contributions: AS*! from the temperature 
dependence of the Bose function ri(a;) = l/(exp(tj/T)— 1) 
and AS2 from the temperature dependence of fermion 
polarization. Numerically we find that the former gives 
a ^ T^ contribution to entropy, while the latter which 
can be written as 



'\<i\<Tuv poo 



AS.= 



_1 r\<i\<J^uv roo o 

—^J d\J dLun{u;)liniD^—UJ,) (11) 



(T^/y^high energy cutoff) gives a singular contribution 
ex — T^lnT. The gauge fluctuation contribution to Xu 
(Ax«) is obtained by taking —d'^/dH^ at H — oi 
AF{H). This approach corresponds to summing the 
bubble diagrams for the vertex correction and the self 
energy correction. It takes the form 



^Xu-- 



1 



iM' 



\ci\<Tuv poo q2 

d^J dujn{uj)Im{D^—^nj^), (12) 



where d^Ii±/dy?p is a short-hand notation for 
9^n_L(w,q;/x_F)/9/i|,|^p=o in which n_L(w,q;^F) is the 
transverse polarization function of Dirac fermions with 
finite chemical potential ^f- This expression, which 
closely resembles that of A 6*2, gives a singular contri- 
bution c>c — TlnT. Note that the expressions for A52 
and Ax„ are also applicable to (nonrelativistic) uRVB 
gauge theory ||l^,|l^ , but they are usually ignored in that 
case because they give only subleading corrections while 
AS*! generates a singular correction c>c T'^/'^ |15|, unlike 
our case in which A52 dominate at low temperatures. 
Summarizing our numerical evaluation, we have 



Ax„ = ^Tln^, 



Act 



vl, 2.6T 



(13) 



at low temperatures (T <'^ Tuv/^) in agreement with 
Eqs. (|). 
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FIG. 1. a)xii of YBa2Cu306.63. The inset: spin Knight 
shifts of YBa2Cu4 0g. The verticallines indicate Tc. Symbols 
are as in Ref. P| and M. Dashed line is the susceptibility Xu 
of free Dirac fermions and solid line is the fit to our thoery, 
which includes gauge fiuctuations. b) "f{T) of YBa2Cua06.67. 
c) Xu of La2-a:Sr3:Cu04 (seetext). d) 7(r) of La2-iSr2,Cu04. 



We now discuss our results in light of the experiments. 
In Fig. la we plot Xu of YBa2Cu306.63j a prototypi- 
cal underdoped (bilayer) cuprate, from the Knight shift 
data of Takigawa et al. M We took the liberty of moving 
the zero of Xu by 0.27 states/eV Cu(2), which is within 
the error bars corresponding to uncertainty in the or- 
bital contributions K'"''' (x„ ex JsT^'P" = K-K'"'''). This 
change avoids the unphysical situation of Ref. Q in which 
63^j"^i7^spm^i7^«pm < Q at T = 0. Further support 

for the adjustment of comes from precision measure- 
ments of the Knight shifts in YBa2Cu408 by Brinkmann 
and collaborators ^ who made substantial upward shift 
of iiT^P*" from their previous values |1^. We find that 
the normal state data of Ref. Q] are well-fitted (solid 
line) by Xu{T;vd,Tuv) = Axu + xl- Here Ax„ is the 
numerical evaluation of Eq. (y_2|) whose low T behavior 
is given by Eq. (O), and Xu is the uniform susceptibility 
of bare Dirac fermions with the same upper cutoff Tuv '■ 

Xu 



= -^TT{Tuvl2T), T[x) = £y/ cosh' ydy. The 
two parameters in the fit are chosen to be vd = 0.76J 
and Tuv = 0.63J, where we set the antiferromagnetic 
exchange energy J=1500 K. We expect the gauge fluc- 
tuations to be suppressed in the superconducting state 
(due to Higgs mechanism) so that Xu should cross-over 
to Xu (dashed line) at low temperatures. This is in qual- 
itative agreement with the data below Tc- The inset of 
Fig la shows a similar fit for the spin Knight shifts of 
YBa2Cu408 P], which is again very good. Thus our the- 
ory can account for the susceptibility in both the normal 
and superconducting states, without the need to adjust 
the energy scale of the gap parameter. Using the same 
parameters V]j and Tuv as in the fitting of YBa2Cu306.63 
Knight shift data, we plot 7 = Acf /T + 7° (where 



7° = ^Q{Tuv/m,Q{x) = !^ vV coBl? ydy) in Fig 
lb. Also shown is the experimental data for 7(T) of 
YBa2Cu30g.67 0- Rough agreement of scales between 
the curves is quite encouraging. 

In monolayer La2_2;Sr2;Cu04, the uniform susceptibil- 
ity is usually deduced from bulk susceptibility by sub- 
tracting the core diagmanetism Xc and Van Vleck para- 
magnetism Xvv Fig. Ic shows Xu of La2-a;Sr2:Cu04 
obtained by subtracting the powder average value Xw + 
Xc — —0.5 states/eV [|l9[ (there's some uncertainty in 
the value of Xw) from the bulk susceptibility x lO- "^^^ 
data can be characterized by Xu = Ax„ -I- x° + Xconst 
with {vD = Qmj,Tuv = 1-17J) for x = 0.10 and 
{vD = 0.79 J,Tuv = 0.65 J) for x = 0.14. Unlike 
the YBCO compounds, temperature independent part 
Xconst > is needed for a reasonable fit. Regarding the 
specific heat data of LSCO, cutoffs significantly smaller 
than the ones used for Xu are needed to fit 7 of the same 
compound in terms of 7 = /S.cf /T + 7" + "/const- In 
Fig. Id we have kept the same vd as in Fig. Ic, but used 
smaller cutoffs {Tuv = 0.8 J for x = 0.1 and Tuv = 0.49 J 
for X = 0.135) to fit the 7-data M. This discrepancy and 
the origin of nonzero jconst and Xconst are not well un- 
derstood. The Xconst > feature in LSCO has been em- 
phasized by some [ pO| to be an important evidence that 
the bilayer structure is important for spin gap behavior. 
Recent experiments on trilayer HgBa2Ca2Cu30g_|_5 |2|] 
and monolayer HgBa2Cu04^^ [ p2| , however, find similar 
spin gap behaviors as in YBCO, suggesting that LSCO 
is a rather special case. 

Despite reasonable agreement, we feel that above com- 
parisons do not provide a conclusive test, because the Tc 
is too high to probe the normal state infrared behavior for 
a wide range of temperature. In fact, most of the bend- 
ing feature seen in the 7(T) data is presumably related 
to the high energy cutoff (the deviation from linear Dirac 
spectrum) which have been treated in a cavalier manner 
here by using a hard cutoff Tjjy. The low-T curvature 
in Xu data {cPxu/dT^ < 0; faster decrease at lower tem- 
perature) seems to support the gauge fluctuation picture, 
but it may not be simple to separate this effect from the 
curvature due to high energy cutoff. Nevertheless, we 
view that the theory advocated here presents a simple 
and appealing picture of the spin gap phase. In this the- 
ory, no new energy scale is introduced to distinguish the 
spin gap phase and the superconducting phase; the Dirac 
velocity in both phases is taken to be the same. Rather, 
it is the gauge fluctuation that distinguishes the phases 
by causing the enhancement of Xu and 7 in the normal 
state. 

Instead of conclusion, we recapitulate some issues that 
have been glossed over. We have ignored the ag field 
whose effect may not be totally innocuous |23]. In fact, 
we have checked that in the absence of bosons the con- 



tributions to Xu and c^' derived from the free energy 
shift due to the ag fluctuation cancel the singular contri- 
butions from transverse gauge fluctuation, in agreement 
with non-renormalization of Dirac velocity in a Lorentz 
invariant situation. Secondly, we have not treated con- 
tributions from the Bose sector, especially in regard to 
the entropy. Lastly, we mention the issue of whether the 
renormalization of fermion propagator feeds back to the 
gauge propagator. In the non-relativistic gauge theory 
|lO,n5[, the density-density correlation function and the 
transverse gauge propagator receive only sub-leading cor- 
rections [1^] . This might not hold any longer in our case. 
At present it is not clear to what extent the transverse 
propagator is modified by "feedback effect" and to what 
extent this affects the physical picture. 
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